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^ : Abstract 

We develop the effective field theory of density fluctu- 
ations for a Newtonian self-gravitating A-body system 
in quasi-equilibrium, apply it to a homogeneous universe 
with small density fluctuations. Keeping the density fluc- 
tuation up to the second order, we obtain the nonlinear 
! field equation of the 2-pt correlation £(r), which contains 

"sj- ; the 3-pt correlation and formal ultra-violet divergences. 

By the Groth-Peebles hierarchical ansatz and the mass 

00 ' 

renormalization, the equation becomes closed with two 



x 



new terms beyond the Gaussian approximation, and their 
coefficients are taken as parameters. The analytic solu- 
tion is obtained in terms of the hypergeometric functions, 
which is checked numerically. With one single set of fixed 
two parameters, the correlation £(r) and the correspond- 
ing power spectrum P(k) match simultaneously the re- 
sults from all the major surveys, such as APM, SDSS, 2df- 
GRS, and REFLEX. The model gives a unifying under- 
standing of several seemingly unrelated features of large 
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scale structure from a field-theoretical perspective. The 
theory is worthy to be extended to study the evolution 
effects in an expanding universe. 
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1 INTRODUCTION 

Great progress has been made in understanding the the large scale 
structure of the universe in past decades. Not only observations 
of the major galaxy surveys such as SDSS [36, 42, 43], 2dF [6, 22, 
26, 29], APM [25, 27], and REFLEX [7, 31] etc, have revealed the 
cosmic structures of increasingly large sky dimension with new de- 
tailed features being found, theoretical studies have also achieved 
important results, through numerical simulations [39, 34], pertur- 
bation method [28, 9, 15, 18, 3, 38], and thermodynamics [30]. 
From view point of dynamics, the Universe filled with galaxies and 
clusters is a many-body self-gravitating system in an asymptotic 
relaxed state, since the cosmic time scale 1/Hq is longer than the 
local crossing time scale [30]. A systematic approach to statisti- 
cal mechanics of many-body systems is to convert the degrees of 
freedom of discrete particles into a continuous field. Thereby, the 
fully-fledged techniques of field theory can apply to study the sys- 
tems [47]. The Landau-Ginzburg theory is a known example in 
this regard. We have formulated such a density field theory of self- 
gravitating systems and applied it to the large scales structure of 
the Universe [44]. Under the Gaussian approximation, the field 
equation of the 2-point correlation function and the solution have 
been derived explicitly. The result qualitatively interprets some 
observational features, but it suffers from insufficient clustering on 
small scales. In this paper, we will go beyond the Gaussian level 
and include nonlinear terms of density fluctuations up to the sec- 
ond order, yielding a more satisfying description of the large scale 
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structure. 

2 NONLINEAR FIELD EQUATION OF COR- 
RELATION FUNCTION 

The universe is represented by a collection either of galaxies, or of 
clusters, including their respective dark halos, as the unit cells with 
random velocities. Although the unit cell, galaxy or cluster, has 
different mass m, both cases correspond to the same mass density 
p(r). We study the asymptotic relaxed state of this Newtonian 
self-gravitating system of N points of mass m with the Hamilto- 
nian H = £f ^ - E? <3 j^=^. Thus the evolution effects will not 
be addressed in this paper. By using the Hubbard-Stratonovich 
transformation [47], the grand partition function of this system 
at temperature T can be cast into the generating functional as a 
path integral Z = J Dfiexp [— (3~ x j d 3 r £(</>)], where is the gravi- 
tational field, (3 = AnGm/c 2 of dimension m _1 , c s = (T/m) 1 ^ 2 the 
sound speed, the effective Lagrangian is C{<f>) = \(xj(j)) 2 — k 2 e^, and 
kj = {AnGpo/ c 2 ) 1 / 2 is the Jeans wavenumber. The term —k 2 e^ has 
a minus sign because gravity is attractive. By the Poisson's equa- 
tion v 2 0( r ) + k 2 e^ = the mass density p is related to the field 
by p(r) = ran(r) = poe^ r \ So po is the constant mass density 
when = 0. We define a dimensionless re-scaled mass density 
field ^(r) = = p(r)/po, and introduce an external source J 
coupling with ip in the effective Lagrangian [44] 

c(:4>,J) = \{^f) — kjip — J 4'. (l) 

J is used to handle the functional derivatives with ease. So far c s 
is the only parameter in place of temperature T, upon which (3 and 
kj depend. The field equation of ip in the presence of J is 

V^-i(v^) 2 + ^ 2 + ^ 2 = 0. (2) 

The n-point connected correlation function is given by G^(ri, . . . , r n ) = 
(^(ri) . . . 5ip(r n )), where 5ip(r) = ip(r) — (ip(r)) is the fluctuation 
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field about the expectation value (ip(r)), and (5i/j) = 0. A stan- 
dard way to derive the field equation of is to take functional 
derivative of the expectation value of Eq.(2) w.r.t. the source J, 
and then set J = [16, 44]. Assuming the large-scale homogeneity 
of the Universe with a constant background density (ip) = ipo, and 
keeping up to the second order of small fluctuation 

# + (^f) , 0) 



i> {i>) + # M V <</>) <</>) 

we obtain the field equation of the 2-pt correlation function 
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V 2 Gf\v) + 2k 2 jG { ? 

1 ._o„,oW_. 1 



^G( 2 )(r)V 2 G( 2 )(0)-^V 2 G?)(0,r,r) 

Wo Wo 

+^-VGf\v) ■ VGf (0)1 = -^ (3) (r), (4) 
Wo J 

where Gf\v - r') = G^(r,r') = (3S(ip(r))j =0 /SJ(r') has been 
used. Eq.(4) is not closed, as it involves the 3-pt correlation func- 
tion G^\ If higher order terms in Sip were allowed in Eq.(3), there 
would be G^\ etc, in Eq.(4). Therefore, we have a hierarchy of field 
equations, typical for the kinetic equation of a generic many-body 
systems. To close Eq.(4), we adopt the Groth-Peebles hierarchical 
ansatz [20] 

(n,r 2 ,r 3 ) = Q[G( 2 Vi,r 2 )Gf (r 2 ,r 3 ) + 
Gf (r 2 , r 3 )G^ (r 3 , n) + Gf (r 3 , n)G^ (n, r 2 )] , (5) 

where the constant Q = 0.5 ~ 1.0 ([15, 11, 23, 24]). Then, Eq.(4) 
becomes closed: 

V 2 Gf } (r) + JfcgG< 2) (r)+a- VG^(r) - 6(VGf(r)) 2 

= -^ (3 V), (6) 
where a = - f )VG( 2 )(0), 6 = Q/ip > 0, and 

^o 2 -2^ 2 + (i,-^)V 2 G( 2 )(0). (7) 

Wo Wo 
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Due to the higher order terms in Eq.(3), the friction term a • VG 
and the nonlinear term b{VG) 2 occur in Eq.(6). As is expected for 
an interacting field theory, 2k 2 is modified by an apparently diver- 
gent term V 2 G^(0). We take k$ as the physical wavenumber like 
in the standard mass renormalization [47]. Eq.(6) is a nonlinear el- 
liptic equation with a point source —^Jq(35^\y). Since 1^(3 oc m/c 2 , 
so galaxies or clusters with greater mass have a higher correlation 
amplitude. This naturally explains why the correlations of clusters 
or of galaxies increase with richness and luminosity [44]. As men- 
tioned earlier, galaxies and clusters are treated on equal footing as 
gravitating particles differing only by their masses, their correla- 
tion functions have the same functional form, differing only in the 
amplitude oc m/c 2 s . This is the observed facts [21, 2]. By isotropy 
of the Universe, one puts Gf\r) = £(r), then Eq.(6) reduces to 

f 'Or) + (- + a)Z'(x) + «*) - b?\x) = -$Pko^, (8) 

where x — k r, £' — and a = - ^)£'(0). Both a and b are 
treated as two parameters. 

3 ANALYTIC AND NUMERICAL SOLUTION 

The Gaussian approximation [5, 44] is recovered by setting a = = 
b and k^ — > 2k 2 in Eq.(6), i.e. keeping only the term ^ ~ in 

Eq.(3). The solution is f(r) oc A c -^^ with A = ^ and k Q = 

( 87T ^ P ° ) 1 , and the power spectrum P(k) = ^ fl | 2 _ , where n is the 

spatial number density. This result qualitatively explains several 
observed features, such as a stronger correlation for more massive 
galaxies, galaxies with a smaller n having a higher P(k) [10, 14] , the 
scaling of "correlation length" tq with the mean cluster separation d 
as ro(d) oc <i°- 3 ~ - 5 [l ? 2, 8, 17, 41], and damped oscillations of £ c (r) 
for clusters with a wave-length 2k jk§ ~ 120 h _1 Mpc [4, 12, 13, 14, 
37, 35]. Here one sees the physical meaning of the sound speed c s . 
Using the background mass density po = pSl m = {S/87iG)HQh 2 Q m 
leads to c s = y^H^htill 2 /k^. Taking the observed periodic length 
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27r//co ~ 100 - 120/r 1 Mpc and htt^ 2 ~ 0.36 by WMAP [32, 33], 
yields c s ~ 1000 ~ 1200km/s, which is roughly that of peculiar 
velocity of clusters or galaxies. Thus, viewing c s as the random 
velocity of clusters or galaxies is qualitatively consistent with the 
observations of the periodicity in £ c (r), of Hq, and of hVt]l 2 . The 
shortcomings of the Gaussian solution are that £(r) is too low at 
small scales, and that P(k) has a sharp peak at k = ko and becomes 
negative for k < ko [44]. 

Now these problems are overcome by the nonlinear Eq.(8), whose 
solution is determined by the boundary condition 

«r c ) = C, £(r c ) = D, (9) 

at some r c . Note that, in fitting with observational data, the am- 
plitude C is higher for clusters than for galaxies, as clarified earlier, 
and the slope D is roughly equal for clusters and for galaxies. The 
range of parameters are taken a = (1.0 ~ 1.3), b = (0.01 ~ 0.05), 
ko = (0.03 ~ 0.06) hMpc -1 . In computation we take r c ~ OA/ko 
h _1 Mpc. Eq.(8) is easily solved numerically. In fact, it has an an- 
alytic solution as follows. By perturbation method, since b <C 1, 
one sets the solution as a series £(x) = E^o^&( x )- Eq.(8) be- 
comes |o + (f + a)ip + £o = -^/^ and & + (| + a)£ + & = g h 
where gi = Ej=o£j&-i-j for i > 0. The solution for £o is a linear 
combination of ([19]) 

V\ = e-^ia, 2; z), y 2 = e~ W(a, 2; z), (10) 

where $ and ^ are the degenerate hypergeometric functions, z = 
(a 2 — 4) 1 / 2 x, and a = 1 + "^Jbj- The solution of ^ is a linear 
combination of Eq.(10) plus a particular solution 

where 2/2) is the Wronskian. The condition in Eq.(9) is satis- 

fied by choosing £ (r c ) = C, £j(r c ) = D, and &(r c ) = g(r c ) = for 
z > 0. As is checked, up to the order z = 4, this analytic solution 
agrees with the numerical one. Once £(r) is known, the power spec- 
trum follows from Fourier transform P(k) = 4tt Jq° £(?") bin / fcr ^ r 2 dr. 
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For galaxies from APM, 2dFGRAS and SDSS, the calculated £, gg (r) 
and P(k) are shown in FIG.l, and in FIG. 2, respectively. For RE- 
FLEX X-ray clusters the calculated £ cc (r) and P(k) are given in 
FIG. 3 and FIG.4, respectively. For SDSS clusters, £ cc (r) is given in 
FIG. 5. It is seen that, for the fixed parameters (a, b) ~ (1.2,0.02), 
and ko = 0.05 h Mpc -1 , the calculated £ gg (r), £cc(r), an d their 
respective P(k) match very well all the major surveys for both 
galaxies and clusters, simultaneously. Thus, our density field the- 
ory gives a decent description of the observational data. 

4 DISCUSSION AND CONCLUSION 

Overall, the solution £(r) of Eq.(8) improves the Gaussian one con- 
siderably, as shown FIG. 6, in which the 2dFGRS active galaxies 
are taken for demonstration purpose. It is found that the nonlin- 
ear term — 6(£') 2 has the effects of strongly enhancing £(r) on small 
scales, and making P(k) flatter in k < ko. The friction term a£' has 
the following effects: slightly increasing the height of £(r) on small 
scales; moving the zeros of £(r) to larger r; strongly damping the 
amplitude of oscillations of £ on large scales [4, 37, 12, 13, 14, 35], 
as seen in FIG. 7; smoothing out the sharp peak of P(k) at ko and 
turning P(k) positive for k < ko. 

The main conclusion of this paper is the following. The universe 
containing galaxies or clusters is viewed as a many-particle system 
in asymptotic relaxed state and can be described by an effective 
density field, whereby the techniques of field theory applies, yield- 
ing a perspective on the large scale structure of the universe other 
than the conventional methods. There appears the Jeans scale 
ko ~ (0.04 ~ 0.06) hMpc -1 , which is the unique scale underlying 
the large scale structure as a gravitational system. Up to nonlin- 
ear terms (5ip) 2 beyond the Gaussian approximation, the nonlinear 
field equation of the 2-pt correlation function of the density fluc- 
tuations has been derived and solved analytically. This analytic 
result of field theory interprets several observed features of large 
scale structures. With fixed values of the parameters (a, 6, ko), the 
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solution matches the observed £(r) and P(k) of both galaxies and 
of clusters, simultaneously. 

Although our results match the observational data on large scales 
very well, our model is still preliminary at the present stage, there 
are several problems that need to be further addressed in the future 
as in the following. 

Firstly the calculated ^ gg (r) for galaxies increases too fast on 
very small scales r < 0.2/i _1 Mpc. This indicates that the model 
may break down on such small scales close to a galaxy size. This 
may suggest that either higher order terms of perturbations should 
be taken into account, or the effects of galaxy formation and local 
virilization need to be included. 

Next, there is a limitation in applying the Groth-Peebles ansatz 
in Eq.(5). As is known, for descriptions of any many-particle dy- 
namic system based on Gibbs-Boltzmann equation, a procedure is 
usually taken, which decomposes the complicated equation into a 
set of differential equations, so that each one in the set is possibly 
manageable. However, thereby, a hierarchy of BBKGY type, or 
the like, is inevitably arises. Different treatments of the hierarchy 
are employed for different systems, and a cut-off of the hierarchy 
is usually is used. For instance, in the case of the photon gas of 
CMB, the multipole decomposition is involved for the temperature 
anisotropics and polarization, and the common practice is to cut 
off the hierarchy of multipoles by letting higher multipole compo- 
nents to zero, yielding a very accurate description of CMB spectra 
[46, 45, 40]). But for the calculation of Gf\r) in our context, one 
can not set G^(ri,r2,rs) to zero, since these are important and 
give rise to nonlinear effects, due to the long range nature of grav- 
ity force. The Groth-Peebles ansatz has been used in our analytic 
treatment, because it is a simple one and an analytic solution can 
be derived. We would like to mention that the ansatz only approx- 
imately reflects the actual distribution since the forms of factor, Q, 
is in fact a function of r. In this case Eq.(6) for Gf\r) would be 
more complicated and analytic solutions for the general case would 
be difficult to derive explicitly. 
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Thirdly, in fitting the observation of galaxies, we have not sepa- 
rated the dark matter from galaxies in our present model. There- 
fore, no bias is introduced and the baryon acoustic oscillations are 
not incorporated in our model. A comprehensive treatment of two 
components, dark matter and baryons, in our theory would require 
substantial extensions to the model discussed the present paper. 

Finally, it should be mentioned that our theoretical model deals 
with only the quasi-relaxed state of the large scale structure of the 
universe, which, as an assumption, is qualitatively good approxima- 
tion since the overall expansion rate Hq is smaller than the particle 
collision rate. It would be much desired to have an extension of the 
present model to take into account of an evolutionary description. 
These issues will to be addressed in our further studies. 
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Figure 1: £, gg (r) fits the galaxy correlation functions of the surveys, APM, SDSS, 
and 2dFGRS, simultaneously. 



Figure 2: P(k) of FIG.l also fits the power spectra of APM, SDSS, and 2dFGRS. 
P(k) oc k' 1 - 5 in (0.05 ~ 0.5)hMp C - 1 . 



Figure 3: £ cc (r) fits the correlation function of REFLEX X-ray clusters ([7]). 



Figure 4: P(k) from FIG. 3 fits the corresponding spectrum of REFLEX X-ray 
clusters ([31]). 



Figure 5: £ cc (r) of SDSS clusters of different richness ([2]) are obtained by the 
same set (a, b) but different initial amplitude. 



Figure 6: The solution £ gg (r) of the nonlinear Eq.(8) improves the Gaussian one 
on small scales and matches 2dFGRAS active galaxies. 



Figure 7: The calculated £ cc (r) with oscillations compared with the observed of 
X-ray clusters ([14, 35]). 
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